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ABSTRACT 

We combine optical and near-infrared AO-assisted integral field observations of the merging 
ULIRG IRAS F17207-0014 from the Wide-Field Spectrograph (WiFeS) and Keck/OSIRIS. 
The optical emission line ratios [N II |/Hn, [S Il]/Ha, and [O I]/Ha reveal a mixing sequence 
of shocks present throughout the galaxy, with the strongest contributions coming from large 
radii (up to 100% at ~5 kpc in some directions), suggesting galactic-scale winds. The near- 
infrared observations, which have approximately 30 times higher spatial resolution, show that 
two sorts of shocks are present in the vicinity of the merging nuclei: low-level shocks dis¬ 
tributed throughout our field-of-view evidenced by an H 2 /Br 7 line ratio of ~0.6—4, and strong 
collimated shocks with a high H 2 /Br 7 line ratio of ~4-8, extending south from the two nu¬ 
clear disks approximately 400 pc (~(y/5). Our data suggest that the diffuse shocks are caused 
by the collision of the interstellar media associated with the two progenitor galaxies and the 
strong shocks trace the base of a collimated outflow coming from the nucleus of one of the 
two disks. 
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1 INTRODUCTION 

Shocks play a significant role in our understanding of galaxy 
evolution. Widespread shocks in galaxy mergers trace galaxy-scale 
gas outflows and inflows that can alter the star formation history 
of a galaxy and their emission can confuse measurements of other 
physical properties in host galaxies. Gas can be shocked at the 
interface between two kinematically distinct gas clouds, wherever 
superso nic motio ns occur. For example, cloud-cloud collisions (e.g. 
ISavedoff. Hovenier & van Leetll 1 9671). tidally-induced gas streams 
jMcDowelletal. 20031: IColina et alJ 2004 : Monreal-Ibero et al.l 


l2006Ll20ld : lFarage et alfe010l : IZakamskall20ld) . and outflows and 

galac tic winds from AGN or star fo r mation (e.g. HeckmanetalJ 


199d: L ehnert & Heckmanl 19961: Veilleu x & Rupke 20021 : 


Lfpari et al.ll2004l : ISharp & Bland-Hawthornll2010h can all produce 


radiat ive shocks to dissipate t heir kinetic energy l Allen et al.l 
120081) . Indeed, the simulations of ICox et al.l (1 20041 . 120061) indicate 
that a significant fraction of the orbital energy in a galaxy merger 
will go into shock-heating the gas, potentially driving additional 
galactic winds. 

With such a wide variety of causes for shocks, naturally 
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their effects can be seen in a number of different types of host 
galaxies, such as galaxy mergers jvan der Werf et al]| 1 9931) . AGN, 
and/or starbursts (e.g. I Veilleux et al.ll2003l ). Shocks can enhance 
a galaxy’s [N Il]/Ha, [S Il]/Ito, and [O I]/Hq emission line ra¬ 
tios by more than 0.5 dex Allen et al.l l2~008f) , which can have 
a number of confounding effects. Comparisons of these line ra¬ 
tios and [O IIl]/H/3, used as diag nostics of the e nergy s ource 
(as in the BPT/VQ87 diagrams: iBaldwin. Phillips & Terlevichl 
1198 it IVeilleux & Osterbrockl 11987h . would misidentify a star- 


Armus et al. 19891; Veilleux et al. 1995; Martin 

1997; Dopita et all 

1997j; Kim et alJ 1998; 

Allen etal. 1999!; 

Veilleux et al. 

1999:; 

Sharp & Bland-Hawthorn 

201 Or Rich et al.1 

201(1 2011, 

2014). 


Measurements that rely on emission line fluxes or ratios, such 
as star formation rate and metallicity, can also be affected unless 
shocks are identified and removed from the calculation. 

Shock-excited molecular gas can also be observed at near- 
infrared wavelengths, with a series of ro-vibrati onal emission 
lines of H 2 in the range ~1.9-2.4 //m dMouril 1 19941) . These 
have been observed both in nearby supernova remnants (e.g. 
[Ohva_, Moorwood & Danziger 1989) an d in extragalactic shocks 
<1 Veilleux et aljfl9 97: Suga iet alll 19971.11999l[2003l: iDavies et aT 


2000; Ve illeux et al] 1 20091: IZakamskal l20ld: Hill_& Zakamska 
































































































2 Medling et al. 


120141 : iDavies et~aD l2014al) . Additional ro-vibrational transitions 
in the mid-infrared have also been used to trace such shocks 
<|Appleton et al.||2 ()0fil : IPeterson et al.ll2012l : iMazzarella et al.ll2012l : 

Iciuver et al .2013). 

Locally, Luminous and Ultra-Luminous InfraRed Galaxies 
(U/LIRGs, with log( L ir /Lq) ^ 11 for LIRGs and ^ 12 for 
ULIRGs) are excellent laboratories to study shocks, as Lir cor- 
relat es with star formation rate, AGN fraction, and merger frac- 
tion (ISanders et a HI 19881: IVeilleux et alJll995l: Zanders & Mirabel! 


1 19961 : Veilleux et al.l 20021 : Ishidal 20041) . each of which can 


U/LIRGs (Runke et al. 20o3; Runke & Veilleux 2005; Rtmke et al. 

2005a 

b 3; Martin 2006; Sturm et al. 2011; Rtmke & Veilleux 

20 lF 

Veilleux et al. 2013l; Teng et al. 20131; Snoon et al. 

2013); 

Rodriguez Zaurin et al. 2013; Arribas et al. 20li Cicone et al. 

20 iF 

and are likely a common cause of shocks. Several stud- 


ies of shocks in nearby galaxi es have linked them t o superwinds 
produ c ed by stellar feedback llHeckman et ail 1990l: iGrimes et all 


I tonkins et al 


2012 


2005; Ric h et al.l I 20 T 1 I: Isot o & Martid l2012l ~fs oto et alj 120121: 


_ _ Hong et al. 201 3l) and AGN-driven outflows 

(Mazzarell a et al.l2012h. Although these sys tems are relatively rare 
in the local universe. iLe Floc’h et alj d20Cd ) have shown that high 
infrared luminosities may have been the norm at redshifts be¬ 
yond 0.7. To understand these cosmically important objects and 
the physical processes that occur in_them, the Great Observatory 
All-Sky LIRG Survey (GOALS: lArmus et alhooF team has com¬ 
piled data across a number of different wavelength re gimes, in¬ 
cludin g Hubble Space Te lescope imaging i n the optical 1 Kim et all 
1201 3l) and near-infrare d lHaan et aljEoiil) . Chandra X-ray obser- 
vations. J lwasawa et al.l l201ll) . Spitzer spectra in the mid-infrared 
dStierwal t et al |2 0 1 3l). GALE X observations in the near and far ul- 
traviolet jHowell et alj 20jC|), and Herschel Sp ace Telescope spec¬ 
tra in the far-infrared l lDIaz-Santos et al.l2013l ). 

Recently, subsamples of the GOALS galaxies have been stud¬ 
ied with integral field spectroscopy, allowing an unprecedented 
look at the physical mechanisms present in U/LIRGs. In lRich et al .1 
120141) . a sample of 27 U/LIRGs was observed with the the Wide 
Field Spectrograph WiFeS to obtain resolved optical spectroscopy 
(R~3000 and 7000 in the blue and red sides, respectively). The 
high spectral resolution afforded by the WiFeS instrument enabled 
the decomposition of emission line ratios and revealed the preva¬ 
lence of sho cks, particularly in late merger stages, dRich et alj20ld 
l201lL 120141) which cause star-forming galaxies to appear as LIN- 
ERs or composite galaxies in a single-aperture analysis. Though 
IFU data are capable of identifying the presence of shocks and their 
relative contribution to the total emission line flux, WiFeS lacks 
the spatial resolution and the long wavelengths necessary to trace 
shocks in the dusty cores of U/LIRGs. 

We present the first results of our new study that combines two 
sets of integral field spectroscopy: optical observations with a wide 
field-of-view and moderately high spectral resolution to identify 
the presence of the shocks, and near-infrared high spatial resolu¬ 
tion observations to investigate the origin of shocks in the cores of 
nearby U/LIRGs. This paper focuses on IRAS F17207-0014, the 
first target for which we have obtained both datasets. 

IRAS F17207-0014 (log (L ir /Lq) = 12.46, 17h23m22.01s 
-00dl7m00.2s, z=0.043, 0.848 kpc arcsec -1 ) is a late-stage ma¬ 
jor merger (merger stage classificatio n 5 out of 6, ‘single or ob¬ 
scure d nucleus with long tidal tails ’: lHaan et al.ll201 ll : [kirn et akl 


|201 3l) . and the most infrared luminous system in the southern sky. 
Although dust obscuration in the optical and limited spatial res¬ 
olution in the near-infrared (-^ ~0 , Z 15) have prevented Hubble 


observations from separating the two nuclei, new adaptive optics 
observations with the OH-Suppressing InfraRed Imaging Spectro¬ 
graph (OSIRIS) hav e shown that the two nuclei are still distinct 
(lMedling~et alj2014l) . though they overlap with a small central sep¬ 
aration in projection (O'/24 ~ 20 0 parsecs). We no te that iMartinl 
ll2006h . lSoto & Martini d2012h . and lSoto et al.l <1 20121) classified this 
galaxy as a double nucleus separated by 2-3 ,, based on kinematics 
measured from their longslit data. Our OSIRIS data reveal that the 
two nuclei are considerably closer, and it is thus possible that the 
second nucleus seen by previous authors was in fact merely a sec¬ 
ond kinematic component (e.g. tidal stream). This galaxy is classi- 
fied as a star forming system based on observations in the optical 
I VeiJleux_emiL[ 1993^ Yuan et all 20KJ), mid-infrared 1 Lutzet_al] 
1 19991 ; Risaliti_et_al. 200fij^StierwidtetalJ 201G; InamietaL 2013 ), 


and X-rav dFranceschini et alj|2003l:llwasawa et alJl201 ll). 

This galaxy has been studied by several other authors to date, 
most finding evidence for outflows on a variety of scales. IMartinl 
d2006l) . using long-slit optical spectroscopy, found Na I D absorp¬ 
tion blueshifted by ~400 km s -1 relative to the Ha emission. In 
a similar analysis, iRupke et alj l2005bl lcl) found Na I D absorp¬ 
tion blueshifted from systemic by ^3 00 km s -1 , also interpreting 
this as outflowing neutral gas. Isturm et~ai1 d201 ll ) looked at Her- 
schel/PACS spectra of the OH 79 /rm line and found evidence of 
blueshifted/outflowing molecular gas with a velocity offset of -100 
km s- 1 (V 85 = -170 km s^ 1 , w m ax = -370 km s -1 ). A number of 
authors have examined optical emission lines for signs of outflows 
as well. lArribas & Colinal d2003l) found no evidence of outflows 
in their INTEGRAL observations (alth ough they do sugges t an in- 
flowing component to the southwest): IWestmoquette et alj I 20 T 2 I) 
also find no evidenc e for outflows in their VIMOS data. How¬ 
ever, (Sotoj&^artinl d2012l) fit multiple components to the emission 
lines in their optical long-slit spectra and produce measurements of 
velocity offset (~-100 km s -1 from systemic) and velocity dis¬ 
persion (~100 km s -1 ) for the shocked components, which they 
sum together to produce an estimated outflow velocity of 200 km 
s -1 ; they also use emission line ratios and s hock models to infer 
shock velocities ranging from 50-600 km s~ x . lArribas et al.ll2014l) 
looked again at INTEGRAL data for this object and found, inte¬ 
grated over the galaxy, an outflowing component to the emission 
lines, with v 0 // set = —151 ± 82 km s -1 and v mal = —306 ± 87 
km s -1 . 

In this paper, we present an updated interpretation of shocks 
in IRAS F17207-0014 from the combination of wide-field optical 
and AO-assisted near-infrared integral field spectroscopy. In Sec¬ 
tion [2] we describe the observations and line-fitting techniques. In 
Section [3] we include the details of our shock models used for 
optical emission line analysis. In Section [4] we present the analy¬ 
sis from both our optical and near-infrared data. In Section [5] we 
discuss the implications of the combined data and in Section [6] 
we summarize the results. Throughout the paper w e adopt Hp = 
70kms _1 Mpc -1 , = 0.28, and Ga = 0.72 dHinshaw et al.1 
■ 20091 ). Physical scales were calculated using Ned Wright’s Cos¬ 
mology CalculatoiQ dWright l [2006l ). 


2 OBSERVATIONS AND LINE FITTING 

As a pilot study for our larger survey, we compare here two sets of 
integral-field observations of IRAS F17207-0014. Our optical inte- 


Available at http : / /www. astro . ucla . edu/ ~wright /CosmoCalc . html 
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Figure 1 . Three-colour image of IRAS F17207-0014 from Hubble ACS 
F435W (blue) ACS F814VF (green) and NICMOS F1&0W (red) from 
ItCim et alj j2013h and lHaan et all<20lH) showing the field of view for the 
WiFeS data (large white box, 25" X 38") and the OSIRIS Kcb data (small 
white box, l'/8 X 2'/7). Images are shown with a square-root flux scale in 
order to emphasize outer structure. 

gral field spectroscopy was taken with the Wide Field Spectrograph 
(WiFeS) as part of a larger survey of shocks in galaxy mergers. A 
more detailed analysis of its shocks and kinematics, along with a 
comparison to the population of mergers, will be presented in Rich 
et al. (in prep). The near-infrared integral field spectroscopy was 
first presented in lMedling et alj d2014h as part of a survey studying 
nuclear kinematics in gas-rich mergers. 

2.1 WiFeS Data 

The first set of data was taken w ith the Wide Field Spectrograph 
(WiFeS: iDopita et al.ll2007l . l2010h on the 2.3-meter telescope at 
Siding Spring Observatory. WiFeS is an image slicer with dual 
beams, enabling red and blue spectra to be obtained simultane¬ 
ously, with different spectral resolutions and wavelength cover¬ 
age. The angular resolution is seeing-limited and is ~1.5" for our 
observations, corresponding to ~1.3 kpc for IRAS F17207-0014. 
The observations were taken on 3 June 2010 UT with the B3000 
and R7000 gratings, giving a spectral resolution of R~3000 in the 
range 3700-5700A and a resolution of R~7000 in the range 5700- 
7000A. These data will be presented in full as part of the WIGS sur¬ 
vey (Rich et al. in prep). The data were sampled with square spatial 
pixels of 0V5 on a side which were then binned 2 x 2 to produce 
1" spaxels. The total field of view is 25" x 38", shown in FigureQ] 
The data were taken in 3 exposures of 1500 seconds each, yielding 
a total on-source time of 75 minutes. The data w ere reduced using 
the WiFeS data reduction pipeline JPonita et alj|2010l). For more 
details on the observations and reduction process, see iRich et al .1 

<201 lh . 

We obtain emission line fits for this study using LZIFU 
(Ho et al. in prep), which_uses the penalized pixel fitting routine 
PPXF dCappellari & Emsellemll2004l) to fit stellar template mod¬ 


els (here we used the libraryof templates at solar metallicity from 
iGonzalez Delgado et al.ll2005l) to the continuum and then fits one, 
two, or three Gaussian profiles to each emission line. LZIFU fits 
both the red and the blue spectra simultaneously, allowing for con¬ 
sistent continuum subtraction across the whole spectrum. By fitting 
both sides simultaneously, LZIFU is also able to fit the same veloc¬ 
ity and velocity dispersion for all emission lines, varying only the 
flux. This fitting method has the beneficial effect of using the high 
spectral resolution of the red side to inform the line profiles of the 
blue side, allowing for high fidelity flux ratios for multiple compo¬ 
nents. Due to signal-to-noise ratio limitations, this dataset was fit 
with the one-component mode only. 


2.2 OSIRIS Data 

The near-infrared set of data was taken on 23 and 24 May 2011 
UT with the OH-Suppressing InfraRed Integral-field Spectrograph 
(OSIRIS; lLarkin et alJl2006D on the Keck II telescope using laser 
guide star adaptive optics (LGS AO). OSIRIS uses a lenslet array 
to obtain up to 3000 spectra at once, at a spectral resolution of 
approximately R~3000. The observations of IRAS F17207-0014 
were taken in two modes: a) using the Kcb filter (1.965-2.381 /rm) 
at the O'/1 spaxeU 1 plate scale, and b) using the Hn4 filter (1.652- 
1.737 /im) at the CT/035 spaxeU 1 plate scale. The field of view of 
the Kcb filter is shown in FigureQ]for comparison; the field of view 
of the Hn4 data is about half that size (see Section l4~2l for relative 
position). Data were taken in ten-minute exposures set in an object- 
sky-object pattern, for total on-source exposure times of 60 and 40 
minutes for the Kcb and Hn4 filters, respectively. 


OSIRIS resolves below the seeing limit by utilizing 


Keck’s laser guide star adaptive optics system (Wizinowich et al. 

2000 

; van Dam et al. 2004; Wizinowich et al. 2006; van Dam et al. 

2006 

). This LGS AO system uses a pulsed laser tuned to the 


Sodium D 2 transition at 589 nm to excite atoms in the sodium 
layer of the atmosphere, around 95 km. The distortions of this laser 
spot are then measured and a deformable mirror is used to flatten 
the wavefront, correcting for atmospheric turbulence in both the 
laser spot and the science data in real time. An on-sky source must 
be available to correct for image motion; this tip-tilt star must be 
brighter than ~18th magnitude in R and within 1’ of the science 
target for the Keck AO system. For these observations, we used a 
nearby star of i?-band magnitude 14.1 and separation of 25V 6. The 
Keck AO system enables diffraction-limited imaging in the near- 
infrared (~(y.'065 in A'-band). 

These OSIRIS data were reduced using the OSIRIS pipeline 
v2.fl using the updated OSIRIS wavelength solution for data 
taken after October 2009 which is now available in version 3 
of the pipeline. For more details on the reduction process, see 
Medli ne et al.l <2014l) . We flux-calibrated o ur Kcb datacube usin g 
the NICMOS A-band image presented in IScoville et alj <200(1) . 
Unfortunately, no equivalent calibration data were available to flux- 
calibrate our Hn4 datacube. 

_ Emissi on lines in the OSIRIS datacubes were fit as in lU et all 

<2013l) and iMedline et al.1 <20 1 4 ). with single-Gaussian compo¬ 
nents. Rather than fit all lines simultaneously, lines were fit accord¬ 
ing to species, allowing for different kinematic structure in ionized 
vs. molecular gas tracers. Here we present morphological results 
from the H 2 emission, which has a series of 5 emission lines in the 


2 Available at http : //irlab .astro.ucla. edu/osir is/pipeline . html 
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A'-band. For more details on the kinematics and other tracers, see 
[Medling et al.l J20141) and U et al. (in prep). 


3 SHOCK MODELS 


In the following sections, we compare our optical emission line 
ratios to shock models, which we describe here. 

The g rid of slow shock models used here was computed by 
Rich et a lTll2010l) fro m the MAPPINGS III code (updated from 


Sutherland & Doniti] 1 1 9931) . Briefly, these models are a nalogous 


to t he fast-sho c k mod els presented in I Allen et al.l 1200 Si) . Follow¬ 
ing Ric h et al.l 1201 lh . we restrict our models to fully-preionized 
shocks with shock velocities ranging from 100 to 200 km s -1 , 
which match our data most consistently. 

The abundance set used in these models is that of 
iGrevesse et al.l d2010h with standard solar dust depletion factors 
from Kimura. Mann & Jessberged d2003l) . We note that using these 
dust depletion factors is valid only for shocks slow enough not to 
sputter the dust grains which may be advected into the shocked 
region. The models were then calculated for a range of metallici- 
ties, ranging from log(0/H)+12 = 7.39 to 9.39; only models with 
log(0/H)+12 = 8.69 or higher were considered for IRAS F17207- 
0014. These models also include a transverse magnetic field con¬ 
sistent with equipartition of the thermal and magnetic fields (B oc 
n,7 * 1 ' /2 ; B = 5fj,G for n c = 10 cm -3 ). We note that thi s choic e 
of magnetic field is consistent with the findings of lThompson et all 
d2006 ). who found that starbursts likely have magnetic fields con¬ 
siderably higher than the minimum magnetic field required to pro¬ 
duce radio emission. The flux ratio of the doublet [S II]6716/[S 
II]6731 is sensitive to magnetic field strength, and thus can be used 
as a check of assumptions. We find that our equipartition choice 
of model predicts a ratio [S II]6716/[S II]6731 of 1.0-1.2, whereas 
shock models using a much smaller magnetic field predict a ratio of 
0.5-0.75. Our observations have [S II]6716/[S II]6731 ~1.3, more 
consistent with a magnetic field in equipartition. 

For shocks of each metallicity and velocity, the emission line 
ratios [O IIl]/H/3, [N Il]/Ha, [S II]/H«, and [O I]/Ha were calcu¬ 
lated for comparison with the observed emission line ratios (see 
Section Ed}. 


4 RESULTS: LINE RATIOS AND MAPS 
4.1 WiFeS Analysis 


The emission line ratios [N II]/H a, [O IIl]/H/3, [S Il]/Ha, and 
[O I]/H« are used as diagnostic tests of the ionization mecha¬ 
nisms present in the gas. Different physical processes, e.g. HII 
regions vs. AGN, will cause different flux ratios of these transi¬ 
tions, and therefore objects fall in different locations in the stan¬ 
dard diagnostic diagrams depending upon their primary source 
of excitation. The hard radiation field from an AGN produces 
more excitations of [OIII], [Nil], [SII] and [OI] compared with 
the thermal ra diation field from the young hot O and B stars 
in HII region s ijBaj dwin et aPll98ll:IVeilleux & Osterbrockjl 19871: 

[ 2 ( 


iKewlev et alJ l200ll : iKauffmann et al] I20031 T Kewlev et al . 20061) . 

Though the diagnostic diagrams are commonly used to distin¬ 
guish between star-forming galaxies and AGN hosts, shocks can 
elevate emission line ratios, causing a spectrum to move from 
the HII region of thes e line diagnostics t owards the compos- 
ite or LINER regions dArmus et all 1 1 9891: IVeilleux et alJ [l995: 


lMartinll997l : lDopita et alJ 19971 : Kim et al.ll998l : lAllen et alJ 199S 


Veilleux et al. 199§ Soto et al. 2012; Hons et al. 20 lit Rich et al. 

2014!!; Newman et al. 20141). Mixing betwee 

n AGN and star for- 
but such mixing in- 
, shifting points to- 
diagrams. By com- 
tially-resolved spec- 
ls, it is thus possible 
y the relative conta¬ 
mination, and AGN 

mation can also exist (Davies et al. 2014bl). 
creases the [O IIl]/H/3 ratio more strongl} 
wards the Seyfert region of the diagnostic 
paring these ionization line ratios in spa 
troscopy to photoionization and shock mode 
to map out the shocked regions and identit 
bution to the luminosity from shocks, star 

(Monreal-Ibero et al. 200Q, 201G; Sharp & Bland-Hawthorn 2010; 

Rich et al. 201Cj, 201 lj; Westmoquette et al. 

201 ll; Freeland et al.! 

20 ll|; Soto & Martin 20121; Davis et al. 2012 

; Fogartv et al] 2012j; 


In Figure [2] we show the standard emission line diagnostic 
diagrams for IRAS F17207-0014 for spaxels where emission lines 
in all relevant species have a signal-to-noise ratio of 3 or higher 
(small black points) and 8 or higher (large black points). We include 
two signal-to-noise ratio cuts to avoid sensitivity to uncertainties 
in noise levels in our spectra. All spaxels fall above the pure star- 
burst classification line in the [N Il]/Ha diagnostic, indicating that 
HII regions alone cannot produce the observed line ratios. We note 
that our line ratios are similar to those measured in the correspond¬ 
ing apertures of the longslit data of ISoto & Martini d2012L their 
Figure 2.28). We have overlaid the shock models (coloured grid 
lines) described in Section [3] Our points form a mixing sequence 
of starburst-dominated (toward the bottom left of the plots and be¬ 
low the Kauffmann pure star formation line) and shock-dominated 
(towards the right of the plots) spaxels. These mixing sequences in- 
dicate a varying f raction of the emission due to starbursts vs. shocks 


cticate a varying traction c 

(Rich et a ll 201(1 [ 20 TTI 1 . 


By combining photoionization models from Starburst99 
jLeitherer et al.HT99f]) and MAPPINGS III with the shock mod¬ 
els discussed previously, we plot four possible mixing sequences 
in Figure [3] varying the metallicities of the models in each case. 
That is, for each metallicity, we plot the range of ionization line 
ratios possible from photoionization models with a range of ion¬ 
ization parameters (log(Q(H)) = 6.5^1e8) combined with shocks of 
a range of velocities (100-200 km s^ 1 ). The lowest and highest 
metallicity models (log(0/H)+12 = 8.69 in the top row and 9.39 
in the bottom row, respectively) do not produce mixing sequences 
able to account for all high signal-to-noise (>8) points, and are 
ruled out. We note that no models are able to reproduce all points 
with moderate signal-to-noise ratios (>3), particularly those in the 
log([0 I]/Ha) plot. This could be due to increased noise in the [O 

I] measurements or to difficulties in the shock models themselves. 
We are in the process of producing new shock models (MAPPINGS 
IV; Sutherland et al. in prep) with updated atomic data that may im¬ 
prove [O I]/Ha ratio predictions. In the meantime, we caution the 
reader against using solely [O I]/Ha measurements as shock diag¬ 
nostics. 

As shown in the second row of panels, our mixing models with 
log(0/H)+12 = 8.99 indicate that the least-shocked spaxels contain 
10% ± 10% shocks, because the lowest spaxels fall between the 
purple and blue lines on each of the diagnostic diagrams. The third 
scenario, with log(0/H)+12 = 9.17, has the lowest spaxels never 
reaching the purple line, and sometimes above even the blue line: 
this suggests that 30% ± 10% of the ionization in these regions 
comes from shocks. It is encouraging that these mixing sequences 
are not highly sensitive to the gas-phase metallicity in this regime. 

We also look at the spatial locations of the spaxels with the 
lowest contribution from shocks in Figure [4] Here we have colour- 
coded the spaxels in the BPT diagrams and then included an image 
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Figure 2. Emission line ratio diagnostic diagrams for IRAS F17207-0014. Likely ionization mechanisms fo r each region of the plo t are suggested based on 
the extreme starburst line (solid r ed line.lKewlev et all200l1 ). the pure star formation line (dashed black line. lKauffmann et alJl2003h . and the Seyfert-LINER 
line (teal dashed line. lKewlev et alI2006h , though these classifications do not take into account the effects of shocks. Line ratios from single-component fits to 
each spaxel with highest sign al-to-noise (>8) are p lotted as large black dots and those of signal-to-noise 3-8 as small dots. Pure shock models from an updated 
version of MAPPINGS III Sutherland & DopitallT9 93; Farage et al. 201 d; iRich et al J|20 1 Ot) are overlaid, demonstrating that shocks likely contribute to the 
line ratios of >80% of these spaxels, causing them to fall above the photoionization line in one or more diagnostic plots. The grid includes shock-induced line 
ratios from shocks with velocities of 100 (yellow, top horizontal grid line), 120. 140, 160, 180, and 200 km s —1 (red, bottom horizontal grid line) and for gas 
with metallicities log(0/H)+12 of 8.69 (blue, vertical grid line on left in first panel and right in others), 8.69, 8.99, 9.17, and 9.39 (green, opposite vertical grid 
line). 


of the galaxy to the left, with each spaxel printed in its matching 
colour and with Ha flux contours to show perspective. (Ha emis¬ 
sion is morphologically similar to the continuum emission, with the 
brightest point indicating the galaxy core.) The colour-coding cor¬ 
responds to shock fraction, with orange spaxels indicating highest 
shock fraction. These diagnostics show that the central regions of 
the galaxy have the highest relative contribution of star formation 
(60-100%), and suggest that shocks are dominant at radii beyond a 
few kiloparsecs, where they contribute 80-100% of the line emis¬ 
sion. Though the relative contribution of shocks (0-40%) is lower 
in the centre, this could be due to an increase in star formation 


transition, 1-0 S( 1) at 2.1218 gm, unless otherwise specified.) 
These vibrational transitions of H 2 can be excited by shock heat¬ 
ing, UV fluorescence or thermal excitation from O and B stars, 
and/or X-ray heating from a nearby AGN l lMour 3 119941 ). Because 
we do not see a broad Br> line, and because n o AGN is evi¬ 
dent from optical 1 Wian_etal. 2010), mid-infrared dStierwalt et al .1 


|201 3l) . or X-ray dlwasawa et al.ll20llh observations, the first two 
scenarios are more likely for this galaxy. If H 2 gas is excited by 
UV light from massive stars, the light from those stars should 
also produce Br 7 emission. Therefore a high H 2 /Br 7 ratio in this 


rather than a decrease in shocks. To test the idea that shocks are 

Goldader et al. 1997; Tecza et al. 

2000). Star-forming regions have 

present in the nucleus, we look to higher spatial resolution data in 

Hs/Brv ratios <0.6 ([Joseph et al 

1984; Moorwood & Oliva 1988‘; 

the central regions. 

Rodriguez-Ardila et al. 2004, 2005; Riffel et al. 2010); shock- 


4.2 OSIRIS Analysis 

To investigate the central spaxels of the WiFeS data in more de¬ 
tail, we turn to the OSIRIS data, with a field-of-view that covers 
only approximately the central 2x2 WiFeS spax els. The GALFIT 
modeling and kinematics presented in lMedling et al.l d2014l) reveal 
that the core of IRAS F17207-0014 contains the nuclei of the two 
progenitor galaxies separated by ~ 200 parsecs in projection. Both 
nuclei appear to host small nuclear disks (R e ff of 200 and 410 
parsecs, respectively). 

Though the A'-band observations are at a lower spatial sam¬ 
pling (O'/1 spaxel -1 ) than the H -band observations, they are par¬ 
ticularly interesting because of the H 2 emission present, which 
traces warm molecular gas. There are five H 2 transitions falling 
within our observed wavelength coverage: 1-0 S(3) at 1.7586 gm, 
1-0 S(2) at 2.0338 /itn, 1-0 S(l) at 2.1218 gm, 1-0 S(0) at 
2.2235 gm, and 2-1 S( 1) at 2.2477 gm. (All five of these transi¬ 
tions are seen in our integrated spectrum, although no individual 
spaxel detects all five transitions. Throughout the rest of the dis¬ 
cussion, when we refer to H 2 flux, we are referring to the strongest 


excited Galactic and extragalactic sources have H 2 /Br 7 ratios 
greater (or much greater) than unity dPuxlev et all 199(3) . The sce¬ 
nario of shock-excited H 2 emission is also supported by key ratios 
of H 2 transitions: 1-0 S(3) / 1-0 S(l) = 0.92, 1-0 S(2) / 1-0 S(0) 
= 1.6, and 2-1 S(l) / 1-0 S(l) = 0.12 all fall within the regime 
predicted by the shoc k models of lShul l & Hollenbach: d 1978h and 
iQiernoff et alJ ((7982) and are similar to those predicted by Mouri 
d 19941) using the model of iBrand et alJ I 1 9891) . We see no differ¬ 
ence in our measured ratios between the high H 2 /Br 7 ratio region 
and the rest of the nuclear regions. 


In Figure[5]we show the emission map of the 2.12 gm H 2 line 
(bottom left panel) and Br 7 line (bottom right panel), with contours 
at higher spatial sampling overlaid. We see that the H 2 emission 
does not trace the Br 7 or continuum emission (top left panel and 
white contours everywhere), and indeed that the H 2 /Br 7 ratio (top 
right panel) confirms that the bulk of the H 2 emission is likely due 
to shocks. The H 2 /Br 7 ratio ranges from ~4-8 below the intersec¬ 
tion of the two disks; we additionally find that most of the field-of- 
view has an H 2 /Br 7 ratio above 0.6, suggesting lower-level shocks 
could be widespread. The H 2 velocity dispersion is slightly but not 
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Figure 3. Emission lin e ratio diagnostic diagram s for IRAS F17207-0014. The extreme starbu rst line (solid red line ,[Kewley et al. 2001), the pure star formation 
line (dashed black line, Kauffmann et al. 2003), and the Seyfert-LINER line (teal dashed line, Kewley et al. 2 0061) are the same as in Figure[2] Line ratios from 
single-component fits to each spaxel with highest signal-to-noise (>8) are plotted as large black dots and those of signal-to-noise 3-8 as small dots. In each 
plot, we show line segments indicating models of 0% shocks (violet) to 100% shocks (brown) in 20% increments, with metallicity log(0/H)+12 = 8.69 (top 
row), 8.99, 9.17, or 9.39 (bottom row). Models with log(0/H)+12 = 8.99 and 9.17 account for all high and most moderate signal-to-noise spaxels. Between 
these two models, the lowest line ratios can be produced by <20% shocks (second row) or 20-^40% shocks (third row). 
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Figure 4. Emission line ratio diagnostic diagram s for IRAS F17207-0014. The extreme starburst l ine (solid red linc.j Kcwley et al. 2001), the pure star formation 
line (dashed black line, Kauffmann et al. 2003), and the Seyfert-LINER line (teal dashed line, Kewley et al. 2( )()(]) are the same as in Figure[2] Line ratios 
from single-component fits to each spaxel with high signal-to-noise (>3 for small points, >8 for large points) are plotted in the left three panels, colour-coded 
according to position in the third plot, with colour divisions evenly spaced in log([0 III]/H/3) + log([0 I]/Ha). On the right, we include a zoomed-in image of 
the galaxy with spaxels colour-coded in the same fashion, and with red contours in Ha flux overlaid to give scale. White pixels indicate regions with signal- 
to-noise < 3 for one or more key diagnostic emission lines. We see that the central spaxels have the strongest contribution from star formation (60-100%), 
with outer spaxels having a higher relative contribution from shocks (80-100%). 


significantly elevated in the shocked region; across the field, most 
spaxels show velocity dispersion >100 km s -1 . 


4.2.1 Characteristics of the Shocked Region 

The region of highest H 2 /Br 7 ratio (~4-8) extends at least 0”4 
(340 pc) from North to South and 170-340 pc from East to West. 

The total H 2 and Bry fluxes from this region are 9.75 x 10 -16 
and 2.28 x 10 -16 erg s -1 cm -2 respectively, producing an aver¬ 
age H 2 /Br 7 ratio of 4.3. These fluxes are not corrected for extinc¬ 
tion. The 1-0 S(l) H 2 line flux we measure is approximately two 
orders of magni tude larger than the larges t flux predicted by the 
shock models of lShull & Hollenbachl dl978l) . IChernoff et alj d 1982h 
found increased H 2 emission produced by shock models including 
a magnetic field; our measured line intensity is about one order of 
magnitude larger than their prediction and covers approximately 
half of their simulated beam. These models were tuned to match 
the shocked molecular hydrogen emission surrounding the Orion 
Molecular Cloud; to our knowledge, a comprehensive grid of shock 
models does not exist. Therefore, the discrepancy in line intensity 
between our observations and these shock models may be due to 
differing shock conditions. The limited number (and, perhaps, lim¬ 
ited applicability to this case) of shock models in the literature pre¬ 
vents us from drawing conclusions about the shock properties from 
this measurement at this time. 

The velocity dispersion of the H 2 emission, integrated over 
the shocked region, is 200 km s -1 ; individual spaxels within the 
shocked region range from 145-200 km s -1 . When integrating 
over the shocked region, the signal-to-noise ratio is sufficient to 
decompose the H 2 profile into two Gaussian components (see Fig¬ 
ure This produces a central Gaussian of v 0 ff se t = —11 ± 25 
km s—1 and a = 159 =t 56 km s -1 containing 36% of the flux 
and a blueshifted Gaussian of v 0 ff se t = —113 ± 31 km s—1 and 
a = 268 d= 24 km s -1 containing the remaining 64% of the flux. 




Figure 5. Top Left: Continuum flux map of the central region of IRAS 
F17207-0014 at high spatial sampling (O'/035 spaxel -1 ) in the HnA 
band (1.652-1.737 /im). The major axes of the two nuclear disks identi¬ 
fied in iMedling et alJ d2014h are overlaid in black for context. Top Right: 
H 2 /Br 7 flux ratio calculated from two bottom panels, with continuum con¬ 
tours overlaid; ratios indicate strong shocks (> 4) to the south surrounded 
by a diffuse background of moderate shocks (> 0.6). Bottom Left: H 2 2.12 
flux map (spatial sampling O'/I spaxel -1 ) with contours from contin¬ 
uum overlaid. Bottom Left: Bottom Right: Br 7 flux map (spatial sampling 
O'/I spaxel -1 ) with continuum contours overlaid. 
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Figure 6. H 2 (black) and Br 7 emission (red) as a function of velocity rela¬ 
tive to systemic (cz=12,900 km s —1 ). The H 2 emission is decomposed into 
two Gaussian fits (dotted purple lines: individual components, dashed pur¬ 
ple line: total fit). Reference velocities are shown for the two disks (from our 
kinematic modeling paper, Medling et al. submitted) along the top and for 
a selection of outflow velocities from the literature along the bottom : traced 
by N a ID absorption (Martin 2006), optical emissi on lines lArribas et al. 
120141) . and OH 79/an absorption iShirm et al.l20nl) . 


5 DISCUSSION 

Optical seeing-limited integral field spectroscopy from WiFeS 
(Section [4. lb . covering the central ~20 x 30 kpc, reveals line ra¬ 
tios indicative of a mixing sequence of shocks. The data are con¬ 
sistent with mixing sequences of metallicities log(0/H)+12=8.99- 
9.17. In both sequences, the central spaxels show the lowest rela¬ 
tive contribution to shocks; the fraction of emission line luminos¬ 
ity contributed by shocks in these cases is 10% ± 10% to 30% ± 
10%, depending on metallicity. We also see evidence of two lev¬ 
els of shocks in the near-infrared AO-assisted integral field spec¬ 
troscopy from OSIRIS (Section l4.2b . indicated by moderate and 
high H 2 /Br 7 ratios (0.6—4 and 4-8, respectively). It is encouraging 
that both sets of data show evidence of shocks in the central region, 
though direct comparison of shock properties is not currently pos¬ 
sible due to the lack of near-infrared shock models. Nonetheless, it 
is clear that shocks penetrate both ionized and warm molecular gas 
in this system. 

Although both datasets show conclusive evidence of shocks, 
it is challenging to derive the shock energetics. The WiFeS data 
presented here lack sufficient signal-to-noise to decompose into 
a shocked and non-shocked components, and we defer the reader 
to I Soto et aP J20 1 2ll for estimates of the energy associated with 
the outflow traced b y opt ical emission lines. In the most extreme 
shocks modeled. iDraine et al .1 dl983h estimates that up to 10% of 
mechanical energy in a shock can be converted to H 2 2.12/im emis¬ 
sion. Using our H 2 luminosity, this puts a lower limit on the me¬ 
chanical energy input rate in the shock of ~3xl0 39 erg s -1 . We 
note, however, that the efficiency of converting mechanical energy 
to H 2 luminosity varies by more than five orders of magnitude in 
their models, based on shock velocity, magnetic fields, and pre¬ 
shock gas density and ionization fraction. Without a comprehensive 


grid of shock models predicting H 2 luminosities and line ratios, we 
cannot constrain this further. 

The complex nuclear dynamics of the central region, resolved 
by our OSIRIS spectra, makes it difficult to determine how the 
H 2 emission relates kinematically to the nuclei. In Figure [6] we 
show the H 2 and Br 7 emission lines summed over the shocked 
region as a function of velocity. We then overplot the nuclear ve¬ 
locities for the East and West disks, from our paper on kinematic 
modeling (Medling et al. submitted), as well as a selection of out¬ 
flow velocities found in previous studies for reference. The range 
of outflow velocities found in the literature spans the H 2 velocities. 
However, it is important to remember that the outflow velocities 
from the literature were found in very different spatial regions of 
the galaxy and/or at poorer spatial resolution; thus, matching ve¬ 
locity profiles is not sufficient for identifying a multiphase outflow. 
We also note that it is not precisely clear what velocity would be ex¬ 
pected from shocked gas resulting from a nuclear collision: it could 
be shocked in place, maintaining a radial velocity between that of 
the two nuclei, or the collisional heating may drive the gas out¬ 
wards, appearing blueshifted (as in the models of ICox et al.ll2004i 
1 20061 ) . 

In the following sections we discuss three possible causes for 
the shocks seen in IRAS F17207-0014: AGN-driven outflows or 
heating, starburst-driven winds, and cloud-cloud collisions from the 
merging of the two progenitor galaxies. We evaluate the evidence 
for each cause for each set of shocks below. 


5.1 Shocks Caused by AGN-Driven Outflows or X-Ray 
Heating 

AGN are capable of driving atomic and molecular gas out 
from the nuclear r egions at moderate (e.g. ~350 km s -1 
in NGC 1266; lA latalo et ahl 1201 ll) to high velocities (~2000 
km s -1 ; ISpoon et al] 201 3U Such outflows are th ought to be 
driv en by radiation pressure from the AGN itself dsilk & Reesl 
Il998h and may be responsible for the existence of black 
hol e scalin g relationsje.g. Sprin ge!. Di Matteo & Hemauisl2005l ; 
iDi Matteo. Springel & Hemquist 120051) . It is therefore plausible 
that shocks in the core of (and, perhaps, throughout) a ULIRG be 
caused by an AGN-driven wind. If an AGN were present, it’s also 
possible that the resulting X-rays could heat the H 2 gas, producing 
the high H 2 /Br 7 ratios that we see. 

In the case of IRAS F17207-0014 we find no evidence for 
an AGN from our data or fr om optical 1 YuanjtiaT 2010), mid- 
infrared dStierwalt et al.ll20lH) . or X-ray ilwasawa et al.ll20lH) ob¬ 
servations. We cannot, however, rule out a Compton-thick AGN ob¬ 
scured to our line-of-sight but affecting gas visible along the plane 
of the sky. If such an AGN were driving an outflow causing these 
shocks we might expect to see high velocity gas as a result - which 
is not present in either our data or in the literatu re ( Martini [20061 ; 

ISoto & Martinll2012l:ISoto et al.ll2012l; iRupke et alj|2005bllch. 

Due to the lack of evidence supporting an AGN in this galaxy, 
we regard this scenario as unlikely. 


5.2 Shocks Caused by Outflows Driven by Star Formation 

Perhaps a more common driver of galactic winds is intense star 
formation activity (see lHeckmanll2003l . for a review). The spaxels 
that are dominated by shocks in our WiFeS data appear to be in the 
outer regions of the galaxy (2-5 kpc) , which is consistent with pre- 
vious observations of galactic winds jRich et all2010l . l201lLl2014l ; 
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ISoto & Marti nl20K3 ; ISoto et al J20 1 3) . The line ratios in the central 
regions are plausibly low not because shocks are not present but be¬ 
cause star formation dominates there (e.g. lArmus et all 1989i ). 

The signal-to-noise ratios of our spectra are generally too low 
to decompose the emission lines into a shocked and non-shocked 
component so it is not possible to say from these data whether the 
shocks we see in the WiFeS data are redshifted or blueshifted rela¬ 
tive to the galaxy, or to measure a shock velocity . How e ver, our re¬ 
sults a re co nsistent with those of [Soto & Martin (20121. [Soto et al.1 
d 2012 h . and lArribas et al J d2014h who also find line ratios increas¬ 
ing with galactic radius in a sample of ULIRGs and attribute them 
to outflows. The additional evidence for outfl ows in both Na I D 
absorp tion l Martin[ 20061 : iRupke et aDI2005bl l3) and OH absorp¬ 
tion dSturm et al.ll201 lli makes it extremely likely that outflows are 
present on large scales. 

The strong shocked region of elevated H 2 /Br 7 flux seen be¬ 
low the two disks could be the origin of such an outflow. As 
this region points south, it is aligned with one of the WiFeS re- 
gi ons of increased shock fraction and consistent with the spectra 
of ISoto & Martini d 2012 ll . which show the strongest shocks to the 
south: this could imply that we are seeing the base of an outflow 
that extends to ~ 5 kiloparsecs. Indeed, the broader H 2 component 
in Figure [t^is blueshifted relative to the East disk with velocity off¬ 
sets similar to the outflow velocities in the literature. Observations 
of other outflows beginning as such collimated structures show that 
they can remain mostly collimated over kiloparsec scales (as in 
iShonbell & Bland-Hawthomlll998lf . Other regions of high shock 
fraction in the WiFeS data could be caused by a second outflow un¬ 
seen in the OSIRIS data, perhaps due to geometry. We calculate the 
mechanical energy rate available in such an outflow as E out fi ow = 
MVout = 5.1 x 10 40 erg s -1 , using M = 4 Mq and v ou t = 200 
km s - 1 from iSoto et al .1 120121) . This available mechanical energy 
rate exceeds the minimum mechanical energy rate (3 x 10 39 erg 
s _1 , see Section^ required to produce our H 2 emission in shocks, 
depending on the shock parameters. 

Most of the OSIRIS field of view shows H 2 /Br 7 ratios of ~2, 
still indicative of shocks, though less strong. In the outflow sce¬ 
nario, one might imagine that uncollimated winds contribute to this 
moderate H 2 /Br 7 level. 

Some galaxies with outflowing winds show a correlation be- 
tween optical emission line ratios an d velocity dispersion (e.g. 
iRich et al.lEoml. 1201 ll: Iho et al Jl2014h . but we see no such trend 
here. It is possible that the lack of a correlation simply means that 
the shock velocity is not related to the turbulent velocity of the gas. 
It is worth noting that the velocity dispersions that we measure do 
not extend down to the few tens of km/s consistent with HH-region 
emission, possibly due to a combination of geometry and/or beam- 
smearing, which may affect our ability to detect such a correlation. 


5.3 Shocks Caused by the Merger of Two Galaxies 

The discovery of two nuclear disks in the core of IRAS F17207- 
0014 so close in projection dMedling et alJlTOUih leads us to con¬ 
sider a third possibility, that some of the shocks seen are due to 
collisions of gas associat ed with the two disks as they are und ergo- 
ing a close passage (as in lHarwit et alJI 1987l : ICox et al.ll2004l) . The 
morphology of the strong (>4) H 2 /Bt '7 region in the OSIRIS data 
is consistent^vith the interface between the two disks. 

IRich et all d201 lh found that the enhanced optical line ratios 
from shocks are easily washed out by star formation, and are thus 
easier to see on the outskirts where star formation rates are lower. 
It is important to note, however, that this does not mean that shocks 


are only present on the outskirts. The situation here is analogous: 
if the disk collision has caused the shocks, we still may not expect 
to see high H 2 /Br 7 ratios directly on top of the two disks, where 
star formation is high and could wash out their signature. Seeing 
a stronger signal of shocks along the interface but below the in¬ 
creased star formation in the disks is consistent with this. 

Given that IRAS F17207-0014 is in a late stage of merger 
(merger class 5 out of 6 : lKim et~al ■ 20l3h. it is plausible to imagine 
that cloud-cloud collisions are not uncommon, and may be widely 
distributed. These may also be responsible for the moderate levels 
of H 2 /Br 7 (~0.6-4) seen in the regions surrounding the nuclei in 
the OSIRIS data, independent of the cause of the highest H 2 /Br 7 ra¬ 
tio. As the OSIRIS field-of-view is completely filled by the overlap 
region of the two progenitor galaxies, such collisions in the inter¬ 
stellar medium would appear as the diffuse background of shocks 
observed. 

The collision of two disks of this size would produce a sub¬ 
stantial amount of mechanical energy. We approximate the energy 
of the system as T = \W = GM di^iM diak2 _ 4 5 x 1Q 58 erg; 
where G is the gravitational constant, Mdiski = 9.8 X 10 9 Mq, 
Mdisk 2 = 4.56 x 10 9 Mq from our dynamical modeling in 
Medling et al. (submitted), and the distance between the two disks 
r=200 pc from lMedline et al.l J2014h : we divide this energy by 100 
Myr as a typical timescale for nuclei to merge, in order to produce 
the average energy input from this source: 1.4 x 10 43 erg s _1 . This 
far exceeds the required amount of energy (3 x 10 39 erg s _1 , see 
Section[5} to produce the measured H 2 emission in shocks, depend¬ 
ing on the shock parameters. 


5.4 Two Mechanisms Contributing to Shocks 

Both the starburst-driven wind scenario (Section \52\ and the 
merger-driven collisional shock scenario (Section 1531 can plau¬ 
sibly describe the data presented in this paper. We propose a hy¬ 
brid scenario, in which the diffuse shocks (H 2 /Br 7 ratio ~0.6^4) 
filling most of our OSIRIS field-of-view are caused by the colli¬ 
sion of the interstellar media associated with the two nuclei, and 
the strong shocks (H 2 /Br 7 ratio ~4-8) to the south of the nuclei 
trace the base of the outflow extending several kiloparsecs, visible 
in the wide-field optical spectroscopy. If the strong shocks do in¬ 
deed show the base of a starburst-driven wind, then deeper optical 
integral field spectroscopy should find the strongest outflow signa¬ 
tures to the south. 

This paper gives a first look at the unique perspective of¬ 
fered by combining a high-resolution view of a galaxy’s core with 
the large-scale information of the system at large. This combina¬ 
tion is an efficient way to look at the fueling and feedback of 
star formation and AGN activity. We are in the early stages of a 
WiFeS+OSIRIS survey of nearby U/LIRGs, enabling us to trace 
spatially-resolved excitation mechanisms from small scales to large 
and study the physical processes involved in star-forming and ac¬ 
tive galaxies and their effects on galaxy evolution. 


6 SUMMARY 

We have presented additional analysis of integral-field spectro¬ 
scopic observations of the ULIRG IRAS F17207-0014, a nearby 
late-stage major merger, making use of wide-field optical emission 
line ratios and complementary high spatial resolution AO-assisted 
near-infrared IFU data. We find: 
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• The optical emission line ratios [N Il]/Ha, [O IIl]/H/3, [S 
Il]/Ha, and [O I]/Ha are elevated relative to pure star formation, 
and show a mixing sequence of star formation and shocks. These 
shocks are wide-spread and have their strongest relative contribu¬ 
tions in the outskirts of the galaxy, which is suggestive of galactic- 
scale winds. 

• High spatial resolution near-IR imaging and kinematics reveal 
that the core of this system is a collision of two nuclear disks. High 
H 2 /Br 7 ratios (4-8) extend from the interface between the disks to 
the south. More moderate H 2 /Br 7 ratios (0.6-4) are present across 
most of the core, indicating the presence of shocks are widespread. 

• We find plausible evidence for two potential causes of the 
shocks seen in IRAS F17207-0014. We propose that the diffuse 
low-level shocks are caused by cloud-cloud collisions driven by the 
merger, while the strong shocks towards the south trace the base of 
the starburst-driven wind seen in both our optical spectroscopy and 
that of other authors. 

• Pairing wide-field IFU data with high spatial resolution IFU 
observations of the central regions is a unique way to characterize 
the physical mechanisms present in the galaxy. We intend to collect 
this combination of data for a large sample of nearby (U)LIRGs in 
order to determine where and how galactic winds are launched, and 
how often they contribute to shocks across the galaxy. 
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